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On 



We give a derivation of the particle-hole RPA equations for 
an interacting multi-fermion system by applying the instanta- 
neous approximation to the amputated two-fermion propaga- 
tor of the system. In relativistic field theory the same approx- 
imation leads from the fermion-antifermion Bethe-Salpeter 
equation to the Salpeter equation. We show that RPA equa- 
tions and Salpeter equation are indeed equivalent. 



equation and RPA equations are indeed equivalent. Con- 



cluding remarks are given in Sec. IV 



II. RPA EQUATIONS AND THE 
INSTANTANEOUS APPROXIMATION 

A. Pole structure of the polarization propagator 
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I. INTRODUCTION 

In the study of systems of many interacting fermions 
one finds that methods used in relativistic field theory 
usually correspond to approximations well-known in non- 
relativistic many-particle theory. An interesting example 
is the correspondence of the Hartree-Fock approximation 
and the so-called Gap-equation used e.g. in the Nambu- 
Jona-Lasinio model jD, which both give a first approxi- 
mation for the many-body problem. An improvement is 
obtained by including correlations in the excited states 
as well as the ground state of the system. In nonrela- 
tivistic many-body theory this leads to the particle-hole 
Random-Phase-Approximation (RPA) equations. 

The RPA equations have been derived by various meth- 
ods (see e.g. ref. The most systematic approach, 
however, is the Green's function method, especially in 
view of possible generalizations. Although the connec- 
tion of the RPA equations to the Bethe-Salpeter equation 
is mentioned in the literature ||, a clear and systematic 
description appears to be missing. It is the purpose of 
this paper to close this gap and to establish the relation 
to the Salpeter equation. 

The paper is organized in the following way: In Sec.|l| 
we will give a derivation of the RPA equations based on 
the Green's function method. The derivation uses two 
approximations, i.e. 

• we apply the instantaneous approximation to the 
amputated two-fermion propagator 

• we substitute the full one-fermion propagators by 
the free ones. 

In relativistic field theory the same approximations lead 
from the fermion-antifermion Bethe-Salpeter equation 
to the Salpeter equation ||. The structure of this 
equation shows many similarities with the RPA equations 
(see e.g. refs. We will show in Sec. [II that Salpeter 



Let H be the hamiltonian describing the dynamics of 
a nonrelativistic or relativistic system of many interact- 
ing identical fermions with ground state \tpo) and excited 
states \ip\), A > 0. The corresponding energies will be 
denoted as E and E\ with E\ < Ey for A < A'. For 
simplicity we assume a discrete energy spectrum, i.e. free 
states are considered in a finite space volume. In the rel- 
ativistic case |^o) is usually the vacuum so that the ex- 
cited states are particle-hole excitations of the vacuum, 
e.g. mesons. 

In order to keep the discussion complete and to intro- 
duce the notation some well-known definitions and facts 
will be recalled in the following. 

Let a a , a) a be fermion field operators with the anti- 
commutator given by {aj,,, a^} = 5 a p. They correspond 
to an orthonormal single particle basis tp a which will be 
specified later. The Heisenberg-picture will be used in 



the following and we define A a (t) 



iHi 



The 



two-fermion propagator in this basis is then given by 

i 2 [G(t,t! ,v,,u')] aa /pp> = 
= (ipo\TA a (t) A a , (f) A\{u) A\, («') |^o) 
= ~{MTA\{u) A a (t) A\,{u') A a , (*') |Vo> (1) 

(compare figJl]). 




FIG. f . The two-fermion propagator G 
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Let u = t + e and u' = t' + e with e > 0. In the limit 

e — > one has 

lim [G(t, t', t + e,t' + e)] aa >pp> = 

= e(t - f) (Vol 4(*) A l>V) l^o) + 
+ e(f - 1) (Vol 4(0 4W l^o) (2) 

and with 1 = J2\ IVaHVaI one obtains 

lim [G(t,t',t + e,f + e)] aa , 0p , = 

IS— *0 

A 

(^oW a a \4>x) (4>x\aj 3 ,a a '\'tpo) + 
+ e(t' - t) e -i(Ex-E )(t'-t) 

(Vol4 a «'l^ A ) (V'Ala^aalVo) (3) 

The term for the ground state (A = 0) can be rewrit- 
ten in terms of one-fermion propagators S^g(t — t') = 

-i(iPo\TA a {t)AUt')\ip Q ) as 



(Vo|o^a Q |Vo) (Vo4' a "'l^o) = 
= -lim^(-e)^,(-e) (4) 

The polarization propagator n(t — t') is now defined as 
lim [G(t,t',t + e,t' + e)} aa > m , = 



- lim S^(-e) S* ,{-e) + i [U(t - t')] aa ,^, 

e— j-0 



(5) 



so that 



i [II(f - t')]aa'W = 
= £ [Q(t_ t >) e -i(^-Eo)(t-t') 



A#0 



(Vola^adVA) (VA|a^aa'lVo) 

+ e(t' - 1) e-^- 5 ") (*'-*) 

(Vo|a^,a a '|VA) (ip\\ala a \ip ) 



With 



0(t)e 



-iBi 



: dv 



1 



2iri I E — v — ie 



(6) 



(7) 



one can write the fourier transform Ii(v) = J dt e %vt IT(i) 
as 

(^o\aJ 3 a a \ip\) (^x\ap,a a i\i/}o) 



E 

A^O L 



v - (E\ - E - ie) 

{i(>o\a' gi a a >\ilJx) (ipx\apa a \4> ) 
v+{E x -E - ie) 



(8) 



Thus the energy levels of the excited states appear as a 
doubled system of poles in the polarization propagator 
IT(^) at v = ±(Ex — Eq — ie) as shown in fig.|| 



Im(v) 
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FIG. 2. Location of the poles of II(^) 

This doubling of the spectrum can be traced back to the 
appearence of the time ordering operator T in the defini- 
tion of the two-fermion propagator G. Any equation that 
looks for the poles of II will therefore obtain a doubled 
eigenvalue spectrum. This statement holds for relativis- 
tic as well as nonrelativistic calculations. It is clear that 
the appearence of the second set of poles does not yield 
any further physical information, since only half of the 
poles can be identified with the eigenvalues of the hamil- 
tonian H . 



B. The instantaneous approximation 

The amputated two-fermion propagator M is defined 
as (see also fig.||) 

[G(t, t', u, u')] aa >/3p> = 

-S^(t - u) S F a ,p,(t' u') + S^,(t - u') S*,p(t> - u) + 
+ EE fdhdhdhdUS^it-t^S^^it'-h)- 

e*i ct2 0.30.4 

■ [M(ti,t 2 ,h,t i )] aia2a3ai S F 3l3 (t^-u)S F il3 ,(t i ~u') (9) 



M 



/ 

2 



FIG. 3. Definition of the amputated two-fermion propaga- 
tor M 

The instantaneous approximation of M used for the in- 
vestigation of one-particle - one- hole propagation is given 
by the ansatz 

M M mst = 6(h - ts) 5(h - U) r(ti - t 2 ) (10) 
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i.e. it is assumed that particles and holes interact instan- 
taneously with each other. 

In nonrelativistic many-body theory this ansatz is 
equivalent to taking into account only diagrams that 
have the appropriate instantaneous structure, as shown 
in fig.|. 



M" 



x- 




FIG. 4. Diagrams that contribute to M inst for a 
4-point-interaction 

The same statement also holds for a relativistic 4-point- 
interaction. In more general relativistic held theories like 
QCD, however, the instantaneous approximation cannot 
easily be interpreted in terms of Feynman diagrams. 

With the relations of the previous section (using t' = 0) 
the instantaneous approximation yields for the polariza- 
tion propagator 

ii& n "{t)] aa ,f l p=SZf l ,(t)SZp{-t)+ 



EE 

Q1Q2 



sL4ti-t)s^,(t 2 ) (11) 



To simplify the notation we define 

Furthermore let 

(AB) a p 1 s 



(12) 



a'3' 



An easy way to represent this definition is to define multi- 
indices as A a p~ ( s ='■ Atary)(5p) =: ^ij s0 that the usual 
matrix multiplication can be applied to Aij . With these 
definitions the Fourier transform of eq. ( pi] ) can be written 

as 



iW™ t {v)=g{v)+g{v)T{v)g{v) 



(14) 



The exact Bethe-Salpeter equation for the amputated 
two-fermion propagator M reads 

M1234 = -^1234 + E E j dt5 dt& dt7 dt& 

Q5Q6 a 7 a$ 

[-^1537 S[ s Sq 7 M 6 284 + K 2 537 5*58 $67 -^6184] (15) 



FIG. 5. The Bethe-Salpeter equation for the amputated 
two-fermion propagator M 

(compare fig.^|) with an appropriate particle-hole kernel 
K and the notation [K(t a , U, t c , td)]a a a b a c a d =■ K abcd 
and analogously for M and S F . 

In a first step one neglects the last (exchange) term. 
The second step is to substitute M mst for M which im- 
plies that also K must have this instantaneous structure, 



K 



K' 



<S(tl-t 3 )*(t2-*4) ^(ti-^) (16) 



After applying the Fourier transformation the instanta- 
neous Bethe-Salpeter equation reads 



T(^) = W{v) + iV{v) g(y) T(u) 



(17) 



Together with i H mst = g + gT g one obtains i H mst = 
g-gV U inst or 

ilP» t {v) = (\g(v)]- 1 -iV{v))- 1 (18) 



C. The RPA-equations 

From eq.(||) we know that IL(v) has poles at v = 
±(E\ — Eq — ie). Onthe other hand we can use eq. (fL8|) to 
obtain a spectral decomposition for W nst {v) in order to 
calculate the pole positions and matrix elements in the 
instantaneous approximation. In the following we there- 
fore look for solutions F(v) of the equation 



([ffMr 1 -iV{v)) F(u) = 



(19) 



(13) To proceed further we approximate the full propagators 
S F in g by the 

free propagators [Sg(t)] a0 = -i ( X o\T A° a (t) [A% (0)]t| Xo ) 



with A° a (t) = e iH ^ a a e- iH °\ where H = £ Q e a ala a 
is some 'free' hamiltonian with ground state |xo)- The 
basis states ip a have thus been chosen as eigenstates of 
Hq, i.e. H(jip a = e a ip a . In the nonrelativistic case an 
appropriate choice for Hq is the Hartree-Fock hamilto- 
nian, whereas in the relativistic case one can use the free 
Dirac hamiltonian with some effective fermion mass to. 

For a system of n fermions, the ground state |xo) of Hq 
is a product wavefunction where the lowest n eigenstates 
are occupied, i.e. a a |xo) = for a > n and a\ |xo) = 
for a < n. In the relativistic case |xo) is the filled Dirac 
sea. With 



Hq al |xo) = (e a + e ) a* |xo) 
H a a |xo) = (-e Q + e ) a a |xo) 



(20) 
(21) 



M 



M 



M 



one finds 



[iS F (t)] a/ 3 = e !£a ' S a p 

■ [Q(t) 9(a - n) - G(-t) Q(n - a)] 



(22) 
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and therefore 



• h is a hermitian matrix and has the structure 



[g° (t)] aa > PP' - [S (*)]o/3' [#0 ( _ *)]a'/3 

■[Q(t)Q(a-n)Q(n- (3) + 
+Q(-t)Q(P-n)Q(n-a)] 



(23) 



In the following we will use capital letters like a = A for 
a > n and small letters like a = a for a < n. Using 
eq.(Q) for the 0-functions and the multi- index notation 
of the previous section we find for the fourier transform 
of 3° 



(s»r r 
wr 1 



(Ab) (B'a 1 ) 



(aB) (b'A') 



-i(eA-£b-v)6( Ab ) {B >a>) (24) 
-i (e B - e a + v) 5 (aB) (b , A , } (25) 



Note that the matrix elements of g are zero if one 
of the multi-indices equals (ab) or (AB). Because of 
iH m3t = g + gTg the same holds for H mst and because of 
ill mst = g — gVU mst also the matrix elements of V with 
these indices are not relevant. It is therefore sufficient 
to consider only the index combinations (Ab) and (bA). 
This simplification is a consequence of approximating the 
full fermion propagators by the free ones. Write 



F = 



(26) 



with Xk = Fr A b) an d Yk — F(bA)- I n the following we 
will assume that X, Y are vectors of dimension D, i.e. 
k = 1 . . . D. Define 2D x 2D -matrices fl and N by 
=:-i(il-vN) so that 

&(Ab)(B'a') = ( e A ~ £b) S(Ab)(B>a>) 
Q(aB)(b>A>) = ( e B - Eo) 5(aB)(b'A') 



N(Ab){B'a>) — b(Ab){B'a>) 
N(aB)(b>A>) = —5(aB){b'A>) 



(27) 



which we will write in a simplified matrix notation as 



_ ( [za- £b) 







(e B -e a ) ) ' J " \ 



■ V = ( J \ I (28) 



Setting ft, := Q + V we find that eq.(|l9|) (multiplied with 
i) can be written as 



hF p 



VpNF p 



(29) 



where p labels the different solutions and F p := F(v p ). 
In the nonrelativistic case one usually approximates the 
kernel V by its lowest order contribution, i.e. the two- 
fermion potential. Then this equation is exactly the 
RPA-equation of e.g. ref. 0. 

Up to now it is not clear how to connect v p and F p with 
the eigenvalues and cigenstates of the full Hamiltonian H. 
It is usefull at this point to recall some of the properties 
of the RPA-equation: 



ft = 



A B 

B* A* 



(30) 



X c 



• Let 

F a = ( ) ( 31 ) 
be a solution with eigenvalue v a . Then 

IT = ( ) (32) 



is a solution with eigenvalue v T = —v a (we as- 
sume that v a is real). For the components of the 
eigenvectors this means that FL^ = [F^^]* and 

r (bA) — l r (Ab)\ ■ 



• If F Pl and F P2 are solutions with v pi ^ v P2 then 

(F P1 \N\F P2 ) := (F P1 ) ] NF P2 = (33) 

• If B is 'small enough' the 2D eigenvalues v p are 
real and nonzero. We will use the index convention 
a = 1 ... D for v„ > and r = D + 1 ... 2D for 
t^y < with v r = —v a . The solutions can then be 
normalized as 



(F p \N\F»)=N p 6 p 



(34) 



with = 1 and iV T = — 1. They form a 2£>- 
dimensional basis, i.e. one can expand a vector i* 1 
as F = Y, p c P F p with c p = N' 1 (F P \N\ F) so that 



i = n; 1 \f») (f"\n\ 



(35) 



The following calculation 

(Nu — ft) -1 N\F P ') = 



(v - N^h)- 1 \F p ) = \F" 

v - Vpt 

E vhr Npl]FP) ( FP \ N \ FP ') 



(36) 



now showes that the spectral decomposition of (Nv — 
ft) -1 is given by 

2D 1 

(Nu hy 1 = E N; 1 \F p ) (F p \ (37) 

Let (a/3) stand for or (aB). Since F ( r a/3) = [i^ a) ]* 
the spectral decomposition for the matrix elements of 
Ii mst (v) can be written as 



4 



M r 



E 

(7=1 



1 



v + v„ 



d 3 p 



1 E 

v 4^ 



(45) 



(38) 



Comparing this equation with the exact pole structure 
given in eq.(|§|) we can identify 



F af3 = (V , 0|4 a «IV' O 

v a = E a — E 



(39) 
(40) 



III. FROM THE SALPETER EQUATION TO THE 
RPA-EQUATIONS 

In relativistic field theory particle-hole excitations of 
the fermionic vacuum are described by the fermion- 
antifermion Bethe-Salpeter equation In the in- 

stantaneous approximation this equation reduces to the 
Salpeter equation From the considerations of the 
previous section we expect the Salpeter equation to be 
equivalent to the RPA equations (p9[). We will show in 
this section that this is indeed the case. 

The Salpeter equation for one fermion flavor can be 
written in the form (see e.g. @,§) 



(41) 



where 



d 3 p' 



-I 



(2tt) 3 
d 3 p' 
(2^)3 



A+(p}[W(p,p')^(p')}A-(p} 
A-(p) [W{p,p') W)\ A + (P) (42) 



with the free Dirac hamiltonian H(p) = 7° (tp* + m), 
the projection operators A ± = (ui ± H{p))/{2uj) and 
uj = \Jp 2 + m 2 (don't confuse H(p) with the full hamil- 
tonian of the previous section). Here M is the mass of 
the bound state, m is the effective fermion mass and W 
is the instantaneous interaction kernel. One can define a 
scalar product by 

(Vi \H = J tr (Vl A+ i> 2 A- - 4 A- V> 2 A+) (43) 

with all quantities depending on p and the notation J = 
J d 3 p/{2irf. 

Let u(p), v(p) be free Dirac spinors (we use the con- 
ventions of H in the following). Define 



w< ri ] (P) ■= U AP) ® v t(~P) 

wilHf) v s {-p) ® ul(p) 



(44) 



We will use the box normalization in the following, i.e. 
we substitute 



Since A + t/j A + — A ip A = we can expand 



r,s=±l/2 

(46) 

with some suitable coefficients b. The factors in front of 
the summation sign have been choosen to simplify the 
notation in the following. 

Solving the Salpeter equation 7i?/> — M tp is equivalent 
to solving 



(47) 



for all given ib\ . 

With cq. (f46|) and the relations for u, v of ref. § we 
compute 



tr 



AVV A+ A- - [«,(+>lt A- wit! A+ 



m - 



tr 



4t¥ A+ w ( Tl A- - \vk->V A" wlTJ, A 



(-)lt A-,„(-) 



9 Orr' Oss' 

m 



tr 

= 
tr 
= 



■,(+)! t A +,„(-) 



A + w y 7> A- - IwiW A" w { 7), A+ 



A+ wit), A- - [ti&Olt A- wit}, A+ 



(48) 



(49) 



(50) 



(51) 



where all quantities depend on p. Using the multiindex 
i = (p, r, s) the scalar product in the box normalization 
can therefore be written as 



1 

-1 



4 +) 



(52) 



We further have 

(il>i\Hih) = bh\T *h) + (ih\V fo) with 
p 

(MVih) = E E tr (55) 



(53) 
(54) 



p p' 



We proceed analogously for the kinetic energy term and 
compute 



5 



tr 



w 



(a] 



Srr' $ss' $aa' 



(56) 



(where a, 



±) so that with u>i — u)(p) 



K +) )*(4 +) ), + (6i 



(4" 




(57) 



For the interaction term we define 



« ■ 1/ u to' 



ir,si 



approximation to the amputated two-fermion propaga- 
tor and by approximating the full fermion propagators 
by the free ones. Our derivation holds for nonrelativistic 
as well as relativistic fermionic systems. Since in rela- 
tivistic field theory the same approximations lead from 
the fermion-antifcrmion Bethc-Salpeter equation to the 
Salpeter equation, this equation should be equivalent to 
the RPA equations. We have shown explicitely that this 
is indeed the case. 

The RPA equations have been analyzed carefully by 
many authors, especially in the context of nuclear physics 
(compare e.g. the references given in g). It would be 
interesting to transfer these results to the Salpeter equa- 
tion, e.g. results on the stability of the RPA or on the 
appearence of spurious solutions. 



with ai,a2 — ± and the multiindices i = (p, 7"i,si), j — 
{$' ,T%, S2). Consider interaction kernels that fulfill the 
relation [W W)V = W(p,p') [ip{p')]' ( (this is 
usually the case for kernels of physical interest). Since 



[w { rV(p)Y 



Wr S (p) we have Vy " 
V i J + = (V^~)* so that we can write 



and 



\,j CLx,Ct2 = =t 



( V++ 



v+- 



&<-> V (v+~y (v++y 




Since are arbitrary the Salpeter equation can now 
be written in matrix form as 



= M 



2uj 
2uj 

1 

-1 



/ v++ v+- 
{(y+-y (v++y 

b (+) 



(60) 
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We identify the eigenvalues of the free hamiltonian with 
the kinetic energies of the free fermions as = — e a = u>i 
and cb = —(b = ^>i- Further we identify the positive 
eigenvalues M with v a = E a — Eq and 



b (+) 
&(") 



F° 



(61) 



Therefore we find that the Salpeter equation (pQ ) has 
exactly the form of the RPA-equations (|2^). 

In this context we would like to mention a work of 
J.Piekarewicz which gives a direct derivation of the 
Salpeter equation using a method similar to our deriva- 
tion of the RPA equations as given in section y. 



IV. CONCLUSION 



In the present paper we have shown that the RPA 
equations can be derived by applying the instantaneous 







